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Abstract 

A  model  of  a  system  including  a  proton  exchange  membrane  (PEM)  fuel  cell  and  its  fuel  processing  section  has  been  developed.  The  goal 
was  to  investigate  the  process  configurations  to  identify  the  optimal  operating  conditions  and  the  optimal  process  structure  of  the  system  by 
applying  modeling  and  process  integration  techniques. 

A  simulation  model  using  commercial  chemical  process  simulation  software  [Belsim-VALI  III,  Version  10:  User  Guide,  Belsim  s.a.,  Rue 
Georges  Berotte  29A,  B-4470  Saint-Georges-sur-Meuse,  Belgium  (2001)]  has  been  developed.  The  model  includes  a  steam  methane 
reforming  (SMR),  a  water  gas  shift  (WGS)  and  a  preferential  oxidation  (PROX)  for  the  fuel  processing  feeding  a  proton  exchange  membrane 
cell  stack  (PEM)  and  a  post-combustion  unit.  The  energy  and  the  material  flows  obtained  by  simulation  have  been  used  to  compute  the  process 
composite  curves  of  the  system.  This  has  been  used  to  identify  the  best  heat  exchange  opportunities  and  to  define  the  optimal  operating 
conditions  of  the  reforming  system  in  order  to  provide  the  best  overall  efficiency  considering  the  balance  of  plant. 

By  improving  the  energy  integration  of  the  system  and  by  using  optimization  of  operating  conditions,  the  efficiency  can  be  raised  from  35% 
for  the  reference  system  to  up  49%  in  the  optimized  design.  This  is  obtained  both  by  process  configuration  modifications  and  by  optimizing 
operating  set-points.  The  model  developed  has  been  used  for  sensitivity  analysis  and  will  serve  as  a  basis  for  further  process  design  including 
the  optimal  structure  determination  and  the  thermo-economic  optimization. 
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1.  Introduction 

Proton  exchange  membrane  (PEM)  fuel  cells  for  station¬ 
ary  applications  are  highly  integrated  systems  including  fuel 
processing,  fuel  cell  itself  and  post-combustion  units.  The 
integration  is  not  only  expressed  in  terms  of  material  flows 
but  also  in  terms  of  energy  balance.  Process  design  has 
therefore  a  big  impact  on  the  performances  of  the  total 
system  [2],  Energy  integration  techniques  have  been  widely 
applied  in  the  chemical  process  industry  to  design  and 
synthesize  new  processes.  The  combined  use  of  process 
simulation  packages  with  energy  integration  techniques 
allows  today  to  design  highly  integrated  systems.  Our  goal 
here  has  been  to  apply  these  methods  to  analyze  and  improve 
the  performances  of  a  total  proton  exchange  membrane  fuel 
cell  (TPEMFC).  As  a  starting  point  for  the  design,  the 
ALSTOM  BALLARD  system  [3]  has  been  used  as  a  reference 
case. 
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The  adopted  methodology  is  the  following: 

(1)  Develop  a  simulation  model  that  allows  computing  the 
energy  requirements  defined  as  a  list  of  hot  and  cold 
streams  as  a  function  of  the  major  degree  of  freedom  of 
the  system,  i.e.  steam  to  carbon  ratio,  reforming 
temperature,  operating  temperature  and  pressure  of 
the  membrane,  fuel  utilization  and  air  excess.  This 
model  may  also  be  used  later  on  to  estimate  equipment 
costs  and  compute  operating  constraints. 

(2)  Compute  the  system  integration  using  well  established 
process  integration  techniques  [4]  into  a  specific 
optimization  frame  work  [5] .  This  step  allows  targeting 
the  integrated  system  performances  without  having  to 
impose  the  heat  exchange  system  design. 

(3)  Optimize  the  system  performance  and  analyze  the 
sensitivity  of  major  decision  parameters  in  order  to 
understand  the  major  integration  bottlenecks.  This  step 
aims  also  at  understanding  the  influence  of  constraints 
and  computing  their  importance  with  respect  to  the 
system  performances. 

(4)  Design  the  heat  exchange  network  for  the  integrated 
system  that  is  able  to  reach  the  targeted  performance. 
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Nomenclature 

'Tcell 

cell  area  (cm  ) 

Fcell 

cell  power  (W) 

Eve 

fuel  cell  power  (W) 

./■add 

7  fuel 

fuel  flow  rate  required  in  the 
post-combustion  subsystem  (kg/s) 

fFC 

7  fuel 

fuel  flow  rate  processed  in  the  fuel 
processing  subsystem  (kg/s) 

fi 

total  molar  flow  rate  of  stream  i 

F 

Faraday  number  (=96485 Cb  moR1) 

hi(T,P) 

molar  total  enthalpy  (including  enthalpy 
of  formation)  of  stream  i  at  a  temperature 

T  and  a  pressure  P 

Cell 

current  density  (A/cm2) 

7cell 

cell  current  (A) 

LHVfuel 

lower  heating  value  of  the  fuel  (kJ/kg) 

n* 

molar  flow  of  component  x  in  stream  i 

N 

the  number  of  the  cells  in  the  stack 

PctW 

cell  pressure  (Pa) 

PROX 

preferential  oxidation  reactor 

S7<+) 

heat  to  be  evaluated  from  (— ),  or  supplied 
to  (+)  unit  u  (W) 

Rcell 

cell  resistance  (Qcm2) 

SMR 

Ti 

steam  methane  reformer 
temperature  of  the  stream  i  (K) 

CNernst 

Nernst  potential  (V) 

^Nernst(^cell) 

Nernst  ideal  voltage  at  Tce\\  (V) 

Tcell 

cell  voltage  (V) 

WGS 

water  gas  shift  reactor 

mole  fraction  of  component  x  in  stream  i 

yx 

1  u 

mean  mole  fraction  of  component  x  in  the 
unit  u:  Yxu  =  (y^+yy/2,  where  i,o 
refer  to  inlet,  respective  outlet  of  unit  u 

Greek  letter 

ni 

conversion  rate  of  component  x  in  unit 

u ■  nl  =  Ku  -  K,u)/nlu 

In  this  paper,  we  will  present  the  methodology  to  realize  the 
first  three  steps  of  the  methodology. 

2.  The  total  proton  exchange  membrane  fuel  cell  model 

The  total  proton  exchange  membrane  fuel  cell  (TPEMFC) 
model  has  been  built  according  to  the  ALSTOM  BALLARD 
system  [3].  It  has  been  divided  in  three  subsystems: 

(1)  The  fuel  processing  (FP)  subsystem. 

(2)  The  proton  exchange  membrane  fuel  cell  (PEMFC) 
subsystem. 

(3)  The  post-combustion  subsystem. 

The  model  has  been  developed  using  a  commercial  modular 
software  (Belsim-VALI  III,  [1])  that  is  using  an  “equation 


Fig.  1 .  Fuel  processing  process  flow  diagram. 


solver”  approach  that  reveals  to  be  very  powerful  and  useful 
for  process  synthesis.  The  calculation  of  the  heat  exchange 
system  using  process  integration  techniques  allows  avoiding 
the  calculation  of  a  complex  integrated  system  with  complex 
heat  exchange  links. 

2.7.  The  fuel  processing  subsystem  model 

The  fuel  processing  subsystem  (Fig.  1)  includes  the 
following  three  main  reactors: 

•  The  steam  methane  reformer  (SMR)  that  proceeds  with 
the  main  reactions  of  steam  reforming  to  produce  hydro¬ 
gen  from  the  methane. 

•  The  water  gas  shift  (WGS)  reactor  that  transforms  the 
remaining  CO  into  additional  hydrogen  at  a  lower  oper¬ 
ating  temperature. 

•  The  preferential  oxidation  reactor  (PROX)  whose  aim  is 
the  conversion  of  the  remaining  CO  into  CO2  to  avoid 
poisoning  the  fuel  cell  membrane. 

The  inlet  mass  flows  presented  in  Fig.  1  are  the  following: 

•  The  stream  #  1  is  the  methane  flow  provided  at  3  bar  and 
25  °C.  Pure  methane  has  been  considered  for  the  calcula¬ 
tion  but  the  model  might  be  easily  adapted  when  real  fuel 
including  other  component  will  be  considered. 

•  The  stream  ^2  is  the  water  flow  required  for  the  steam 
methane  reformer  whose  flowrate  is  controlled  by  the 
steam  to  carbon  ratio. 

•  The  stream  ^6  is  a  pure  oxygen  flow  required  for  the 
PROX. 

As  this  model  is  developed  for  system  design,  the  pressure 
drops  have  been  neglected  and  the  operating  pressure  has 
been  fixed  at  3  bar. 

2.1.1.  Steam  methane  reformer 

The  steam  methane  reformer  performs  the  following 
reactions: 

•  The  endothermic  steam  methane  reforming: 

CH4  +  H20  ^  3H2  +  CO,  A77r  =  206. 11  kJ/mol. 

•  The  exothermic  water  gas  shift: 

C0  +  H20  ^  C02 +  H2,  A77r  =  — 41.2kJ/mol. 
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The  two  reactions  are  supposed  to  be  near  equilibrium. 
A  temperature  difference  with  respect  to  the  equilibrium  is 
used  to  represent  the  imperfection  of  the  catalyst.  We 
assume  that  the  size  of  the  catalyst  bed  will  be  sufficient 
to  approach  equilibrium.  The  reactor  is  supposed  to  be 
isothermal,  meaning  that  heat  has  to  be  supplied  to  maintain 
the  temperature  in  the  reactor  from  an  external  energy 
source.  In  order  to  define  appropriately  the  energy  require¬ 
ment  of  the  heat  exchange,  we  have  assumed  that  the  inlet 
streams  will  be  preheated  up  to  the  reaction  temperature. 
This  assumption  is  indeed  conservative  in  the  sense  that  it 
does  not  take  account  for  the  reactor  temperature  profile  that 
will  be  below  the  temperature  profile  obtained. 

In  this  study,  the  temperature  of  the  reformer  T0  smr 
and  the  steam  to  carbon  ratio  have  been  considered  as 
being  the  major  decision  variables  for  the  optimization. 
Other  parameters  might  be  also  chosen  like  the  pressure 
or  the  size  of  the  catalyst  beds  for  further  thermo-economic 
optimizations. 


the  reactor  by  heat  exchange  (see  the  process  integration 
chapter  for  comments  concerning  the  approach). 

2.2.  Proton  exchange  membrane  fuel  cell  model 


2.1.2.  Water  gas  shift  reactor 

We  consider  a  WGS  system  configuration  with  only  one 
water  gas  shift  reactor  that  proceeds  with  the  exothermal 
WGS  reaction  (CO  +  HoO  ^  COo  +  H2)  that  is  supposed 
to  be  near  the  equilibrium.  In  order  to  proceed  with  only  one 
reactor  the  temperature  of  the  system  as  been  fixed  at  250  °C 
using  new  catalyst  as  mentioned  in  [6].  The  isothermal 
conditions  will  be  maintained  by  heat  exchange  in  the 
reactor.  This  means  a  special  design  for  the  reactor  that 
will  differ  from  the  conventional  adiabatic  configurations. 
The  isothermal  conditions  are  assumed  in  order  to  avoid 
energy  losses  in  the  heat  source  definition  and  guarantee  the 
lower  possible  equilibrium  temperature.  This  approach  is 
similar  to  the  approach  used  in  the  reformer  but  for  an 
exothermal  reaction.  The  reactor  inlet  is  supposed  to  be 
cooled  down  to  the  reactor  temperature  before  entering  the 
reactor.  This  defines  a  hot  stream  that  will  be  considered  in 
the  energy  integration  calculation. 

2.1.3.  Preferential  oxidation  reactor 

The  preferential  oxidation  reactor  is  used  to  eliminate  the 
CO  that  has  not  been  converted  in  the  WGS  reactor.  This 
reactor  is  required  to  reach  very  low  level  of  CO  content  in 
the  fuel  stream  at  the  fuel  cell  inlet  in  order  to  avoid  poisoning 
of  the  membrane.  The  preferential  oxidation  is  one  of  the 
possible  technology  for  this  task.  It  is  using  pure  oxygen  to 
proceed  with  the  following  reaction:  CO  +  02^  CO2. 
Unfortunately,  the  selectivity  of  the  catalyst  will  not  avoid 
the  combustion  of  some  hydrogen  in  the  gas  stream  with  the 
following  reaction:  H2  +  (1/2)02  ^  IUO.  In  the  model,  the 
CO  content  at  the  reactor  outlet  is  fixed  at  10  ppm.  The  pure 
oxygen  flow  rate  (stream  #10)  is  computed  as  a  function  of 
the  CO  flow  rate  assuming  one  mole  of  O2  per  mole  of  CO 
according  to  the  results  published  by  [7],  An  isothermal 
operation  at  a  temperature  of  200  °C  has  been  considered 
for  the  PROX  reactor,  if  needed,  heat  will  be  removed  from 


The  PEMFC  subsystem  has  been  built  considering  sepa¬ 
rately  the  cathode  and  anode  sides.  A  schematic  of  the 
PEMFC  submodel  is  given  in  Fig.  2. 

The  operating  pressure  of  the  fuel  cell  is  3  bar.  In  the 
integrated  system,  the  air  compression  work  will  be  pro¬ 
duced  by  expanding  the  post  combustion  gases  in  an  expan¬ 
der.  The  compressor  isentropic  efficiency  (iff)  is  86%.  The 
inlet  air  condition  are  1  bar,  25  °C  and  a  relative  humidity 
of  60%. 

The  equations  of  the  PEMFC  are  the  following: 

•  Nemst  potential  (V): 


UNernst  —  U 


0 

Nernst 


(rcell)+^ln| 


yH2 

1  anode 


(if2 


d/2' 


anode  > 


IF 


/H2o 

anode 


RFcell  (  ^ 

+  ln  cell) 


•  Cell  voltage  (V).  Weil  —  fAjcmst  /Aeidccii  • 

•  Hydrogen  flow  rate  through  the  membrane:  ftfc2athode  =  'ff- 

•  Cell  current  (A):  7ceii  =  Ceil  Weil- 

•  Cell  power  (W):  Pceii  =  Weil  Ceil- 

•  Fuel  cell  power(W):  Ppc  =  WeiW 

•  Hydrogen  oxidation  at  the  cathode:  H2  +  (1/2)02  H2O. 

•  Energy  balance: 


inlets  outlets 

fjhj{ Tcell ,  P cell )  —  foK( Well)  Well)  =  EpC  +  Q^C 

i  =  1  0=  1 

•  Fuel  utilization:  r/comb  =  (uf/node  -  n^2anode)/nfa2node  The 
fuel  utilization  in  the  cell  (t]com b)  has  been  set  at  95% 
according  to  [8].  This  means  that  5%  of  the  hydrogen  at 
the  inlet  of  the  cell  is  not  used  in  the  fuel  cell  and  will  be 
used  as  depleted  fuel  to  balance  the  heat  requirement  of 
the  fuel  processing  section.  The  sensitivity  of  this  para¬ 
meter  will  be  described  later  on. 
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The  PEMFC  subsystem  is  supposed  to  be  isothermal  and 
isobar.  The  cell  temperature  Tceu  will  be  considered  as  a 
decision  parameter  for  the  optimization.  Heat  exchange  is 
realized  to  maintain  the  cell  temperature.  We  assumed  that 
the  hydrogen  combustion  is  complete  on  the  cathode  side 
and  that  only  the  hydrogen  protons  are  crossing  the  mem¬ 
brane  [8].  Water  migration  through  the  membrane  has  been 
ignored  since  it  might  go  in  both  direction.  Water  manage¬ 
ment  strategies  have  not  been  considered  in  the  model.  The 
sensitivity  analysis  procedure  will  be  used  to  evaluate  the 
possible  water  management  strategies  and  their  impact  on 
the  system  efficiency. 

The  stream  at  the  cathode  outlet  is  cooled  down  and  the 
water  is  separated  to  be  partly  recycled  in  the  steam  methane 
reformer. 

2.3.  Post-combustion  subsystem  model 

The  depleted  fuel  of  the  PEMFC  subsystem  is  burned  in 
the  post  combustion  subsystem  (Fig.  3).  The  generated  heat 
will  be  used  to  balance  the  energy  requirement  of  the  fuel 
processing  section.  Supplementary  firing  of  natural  gas  will 
be  considered  if  the  energy  content  of  the  depleted  fuel  is  not 
sufficient  to  close  the  balance.  The  supplementary  fuel 
requirement  will  be  computed  using  the  process  integration 
techniques  as  described  in  the  next  section. 

The  inlet  streams  #8  and  #12  are,  respectively,  the 
cathode  and  anode  outlets.  The  assumptions  are  the  follow¬ 
ing: 

•  The  combustion  is  complete. 

•  Stoechiometric  combustion  has  been  assumed.  This  is  not 
the  case  in  reality,  this  assumption  is  made  to  compute  the 
real  potential  in  terms  of  energy  integration.  In  reality,  the 
air  excess  is  related  to  the  water  management  strategy:  an 
increase  of  the  air  flow  rate  allows  to  reduce  the  liquid 
water  production  on  the  cathode  side,  but  in  the  same  time 
it  reduces  the  post-combustion  temperature  decreasing 
by  the  way  the  available  heat  for  the  reforming.  When 
designing  the  heat  exchange  system,  the  sensitivity  of  the 
assumption  will  be  tested  with  respect  to  the  system 
efficiency  and  might  be  removed. 

•  The  depleted  fuel  and  air  preheating  temperatures  will  be 
considered  as  an  optimization  parameter  for  the  process 
integration. 

•  The  flue  gases  will  exchange  heat  with  the  other  sub¬ 
systems  before  being  expanded  in  a  turbine  to  produce  the 


Fig.  3.  Post-combustion  subsystem  process  flow  diagram. 


air  compression  work.  The  inlet  temperature  in  the  turbine 
will  be  computed  to  balance  the  air  compression  work 
requirement. 

•  After  expansion,  a  new  heat  exchange  is  performed  to 
cool  down  the  streams  to  a  stack  temperature  of  120  °C 
and  a  pressure  of  1  bar. 

2.4.  Heat  exchange  subsystem  model  using  process 
integration  techniques 

Intentionally,  the  heat  exchangers  have  not  been  consid¬ 
ered  in  the  simulation  model  described  above  because 
the  heat  exchange  system  has  been  considered  as  being 
unknown  apriori.  Our  goal  is  to  define  the  heat  exchange 
network  structure  for  optimized  operating  conditions.  We 
have  therefore  used  the  hot  and  cold  composite  curves 
concept  of  the  pinch  technology  to  model  the  integrated 
heat  exchange  system  without  having  to  impose  a  heat 
exchange  network  structure.  In  this  approach,  the  results 
of  the  simulation  model  will  characterize  the  hot  and  cold 
streams  of  the  system. 

The  process  integration  method  (or  pinch  technology) 
(e.g.  [4])  has  been  developed  for  analyzing  the  energy 
integration  in  industrial  processes  to  help  identify  energy 
savings  opportunities.  This  method  is  typically  applied  in 
two  major  steps.  The  first  is  a  targeting  phase  that  aims  at 
identifying  the  minimum  energy  requirement  of  the  system 
by  computing  the  possible  energy  recovery  from  the  hot 
streams  to  heat  up  the  cold  streams  of  the  system.  The 
second  step  is  the  synthesis  of  the  heat  exchange  network 
that  will  implement  the  targeted  energy  recovery.  This 
method  has  been  widely  applied  in  the  process  industry 
for  the  both  grassroots  and  retrofitting  problems.  It  has  also 
been  widely  used  as  a  tool  for  process  synthesis.  Process 
integration  studies  start  with  the  definition  of  a  list  of  hot  and 
cold  streams.  The  hot  and  cold  streams  define,  respectively, 
heat  sources  and  heat  requirements  of  the  system  that  are 
characterized  by  a  heat-temperature  diagram  usually  defined 
by  a  heat  load,  an  inlet  and  a  target  temperature.  The  heat 
sources  are  then  composed  to  compute  the  hot  composite 
curve  that  represents  the  heat  availability  in  the  system  as  a 
function  of  the  temperature.  The  same  procedure  is  applied 
for  the  cold  streams  to  draw  the  cold  composite  curve.  The 
composite  curves  can  be  considered  as  being  one  single  hot 
and  cold  streams  that  will  exchange  heat  in  a  virtual  counter- 
current  heat  exchanger  that  models  the  heat  exchanger  net¬ 
work  of  the  system  under  study.  Considering  that  the  heat 
exchange  will  be  technically  feasible  if  the  temperature 
difference  between  the  hot  and  cold  streams  will  always 
be  superior  to  a  predefine  A7jnm,  the  maximum  heat  recovery 
by  heat  exchange  between  the  hot  and  the  cold  streams  will 
be  obtained  when  the  A7jnm  constraint  will  be  activated.  This 
point  is  called  the  system  pinch  point.  By  heat  balance,  one 
may  then  compute  the  minimum  energy  requirement  of  the 
system  and  the  minimum  heat  losses  to  be  evacuated  from  the 
system.  When  some  of  the  streams  have  unknown  flow  rate  or 
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Table  1 


The  Fuel  processing  subsystem  streams,  syngas  composition:  CO, H20, H2,  CO2,  CH4 


# 

Stream 

Description 

^target 

Flow  type 

1 

Cold 

Water  preheating  for  SMR  reaction  (including  liquid  water 
preheating,  the  vaporization  and  the  steam  superheating) 

20  °C 

Z.SMR 

h2o 

2 

Cold 

Methane  preheating  for  SMR  reaction 

25  °C 

Z.SMR 

ch4 

3 

Hot 

The  stream  between  the  SMR  and  WGS  reactions 

Z.SMR 

Ti,  WGS 

syngas 

4 

Hot 

The  stream  between  the  WGS  and  PROX  reactions 

To,  WGS 

Z.PROX 

syngas 

5 

Hot 

The  generated  heat  by  the  WGS  reaction 

T0,  WGS 

To, WGS 

Heat 

6 

Hot 

The  generated  heat  by  the  PROX  reaction 

Z.PROX 

To,  PROX 

Heat 

7 

Cold 

The  required  heat  by  the  SMR  reaction 

Zi.SMR 

To, SMR 

Heat 

Table  2 

The  PEMFC  subsystem  streams,  flue  gas  composition:  H20,  C02, 02,  N2 

# 

Stream 

Description 

T 

1  in 

T 

1  target 

Flow  type 

8 

Hot 

The  stream  between  the  fuel  processing  subsystem  outlet  and  PEMFC  inlet 

A, PROX 

Zell 

Syngas 

9 

Hot 

The  air  compressor  outlet  cool  down  to  the  cathode  inlet 

To,  comp 

Z,cell 

Air 

10 

Hot 

The  stream  between  the  cathode  outlet  and  the  water  separator 

To,  cell 

35  °C 

Flue  gas 

11 

Hot 

The  generated  heat  by  the  PEMFC 

Zell 

Zell 

Heat 

Table  3 

The  post-combustion  subsystem  streams 

# 

Stream 

Description 

Zn 

T 

1  target 

Flow  type 

12 

Cold 

The  air  preheating  for  the  combustion 

35  °C 

500  °C 

Flue  gas 

13 

Cold 

The  depleted  fuel  preheating  for  the  combustion 

Zell 

500  °C 

Syngas 

14 

Hot 

The  generated  heat  by  the  combustion 

To,  comb 

Z.exp 

Flue  gas 

15 

Hot 

The  combustion  gases  stream  after  the  expander 

To, ex  p 

120  °C 

Flue  gas 

temperature,  optimization  methods  are  used  to  compute  the 
optimal  values  that  give  the  best  integration  (e.g.  [5]). 

Using  the  composite  curve  calculation  together  with  a 
simulation  model,  we  define  the  simulation  of  the  heat 
exchanges  network  without  knowing  about  its  configuration 
and  we  will  be  able  to  compute  the  influence  of  the  major 
decision  variables  to  later  on  define  the  best  system  config¬ 
uration.  The  difficulty  here  is  that  the  heat  exchange  model 
has  to  be  solved  as  an  optimization  problem. 

The  efficiency  of  the  application  of  the  pinch  technology 
is  highly  related  to  the  appropriate  definition  of  the  hot  and 
cold  streams,  i.e.  the  definition  of  the  heat  load  and  the 
temperatures  in  the  system.  These  have  to  be  defined  in  such 
a  way  that  heat  recovery  options  will  not  be  missed  or 
hidden.  This  explains  the  choice  of  isothermal  reactions  and 
some  of  the  preheating  operations  that  have  been  introduced 
in  the  model.  Once  the  results  will  be  obtained,  the  useful¬ 
ness  of  the  heat  exchange  requirements  will  be  analyzed  and 
only  the  useful  exchanges  will  be  considered  in  the  heat 
exchanger  network  design. 

Tables  1-3  show  the  definition  of  the  hot  and  cold  streams 
of  the  system. 

The  heat  sources  and  sinks  are  used  to  compute  the  hot 
and  cold  composite  curves  (Fig.  4).  Assuming  a  AZmin  of 
10  °C  and  25  °C  at  the  SMR,  we  have  computed  the  heat 


recovery  between  the  heat  sources  (hot  streams)  and  the  heat 
requirements  (cold  streams). 

The  A  Zinin  constraint  is  used  to  impose  the  technical 
feasibility  of  the  heat  exchanges.  The  pinch  point  is  the  tem¬ 
perature  at  which  the  constraint  is  activated.  If  the  available 
heat  in  the  hot  streams  above  the  pinch  point  is  lower  than 
the  heat  required  by  the  cold  streams,  a  supplement  of  heat 


Combustion 


Fig.  4.  The  composite  curves  of  the  TPEMFC  system  with  a  steam  to 
carbon  ratio  of  1.4  (SMR  and  PEMFC  temperatures  fixed  at  750  and  85  °C 
and  a  fuel  utilization  of  95%)  and  a  stoechiometric  combustion. 
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will  have  to  be  supplied  to  the  system.  In  our  case,  the  heat 
supplement  will  be  supplied  by  a  supplementary  firing  with 
natural  gas.  The  combustion  will  define  a  hot  stream  (flue  gas 
from  the  adiabatic  temperature  of  combustion  to  the  stack 
temperature).  It  should  be  noted  that  if  the  pinch  point 
temperature  is  the  reforming  temperature  (high  temperature), 
the  supplementary  firing  should  take  into  account  the  pos¬ 
sible  air  preheating  below  the  pinch  point  that  will  reduce  the 
fuel  consumption.  In  the  model  developed,  we  used  an 
optimization  model  [5]  that  allows  the  calculation  of  the 
optimal  fuel  flowrate  and  the  corresponding  optimal  air 
preheating  temperature. 


3.  Results  and  discussions 

The  integrated  process  model  has  been  used  to  determine 
the  optimal  operating  conditions  to  be  considered  in  the 
system. 

The  influence  of  the  following  parameters  has  been 
studied: 

•  The  steam  to  carbon  ratio  (steam/carbon). 

•  The  SMR  temperature  (T0,smr)- 

•  The  PEMFC  temperature  (rcen). 

•  The  fuel  utilization  (rjcomh). 

Analyzing  the  composite  curves  of  the  system  allows  then  to 
study  the  process  modifications  that  will  further  improve  the 
overall  efficiency  of  the  system. 

3.1.  The  self-sufficient  limit 

One  of  the  interests  of  the  combined  use  of  process 
simulation  and  process  integration  techniques  is  the  deter¬ 
mination  of  the  self-sufficient  limit  of  the  TPEMFC  system 
that  is  defined  as  being  the  limit  at  which  the  depleted  fuel 
at  the  outlet  of  the  fuel  cell  allows  satisfying  the  energy 
requirement  of  the  system.  This  limit  is  obviously  related 
to  the  decision  parameters  that  have  been  used  in  this 
study.  Figs.  4  and  5  show  the  composite  curves  obtained 
for  two  different  values  of  the  steam  to  carbon  ratio  at  fixed 
T0, smr  =  750  °C,  Tcei I  =  85  °C  and  t]comh  =  95%.  In  Fig.  4, 
the  depleted  fuel  satisfies  the  energy  requirement  of  the 
system.  When  increasing  the  steam  to  carbon  ratio,  a  pinch 
point  is  created  by  the  steam  production  (at  the  vaporization 
temperature),  because  the  heat  requirement  has  been 
increased  and  because  the  excess  of  water  leads  to  a  higher 
water  content  in  the  depleted  fuel  and  also  because  the  air 
flowrate  reduces  the  adiabatic  temperature  of  combustion.  In 
Fig.  5,  we  should  also  observe  an  increase  of  the  WGS 
reactor  performance  that  corresponds  to  a  decrease  of  the 
PROX  reactor,  revealing  that  the  CO  conversion  has  been 
higher  in  the  WGS  reactor  due  to  the  presence  of  an  excess 
of  water  in  the  syngas  produced  by  the  reformer.  We 
should  also  remark  that  by  the  non  selectivity  of  the  PROX 
catalyst,  part  of  the  oxygen  entering  the  reactor  will  burn 


Fig.  5.  The  composite  curves  of  the  TPEMFC  with  a  steam  to  carbon  ratio 
of  2.4  (SMR  and  PEMFC  temperatures  fixed  at  750  and  85  °C  and  a  fuel 
utilization  of  95%)  and  a  stoechiometric  combustion. 


part  of  the  hydrogen  produced  reducing  the  production 
of  electricity  in  the  cell.  The  results  of  Fig.  4  are  obtained 
for  a  stoechiometric  combustion  of  the  depleted  fuel.  If  air 
excess  is  used  to  control  the  water  balance  of  the  cathode, 
we  obtain  the  composite  curves  of  Fig.  6.  The  combustion 
temperature  is  reduced  and  the  pinch  point  is  activated  at  the 
reformer  temperature.  In  the  case  of  Fig.  6,  the  air  excess 
correspond  to  a  near  saturation  condition  at  the  cathode 
outlet. 

When  the  system  is  not  self-sufficient,  additional  natural 
gas  is  used  to  balance  the  system  requirement.  In  this 
case,  the  overall  efficiency  of  the  system  is  computed  by 
Eq.  (1). 


LH  Vtel(/«  +fS) 

where  Pc\cc  is  a  generated  electric  power  (kW),  FHVfuei  the 
lower  heating  value  of  fuel  (kJ/kg),  the  fuel  flow  rate 
processed  in  the  fuel  processing  subsystem  (kg/s)  and/^)1  is 
the  additional  flow  rate  required  in  the  post-combustion 
subsystem  (kg/s). 

Fig.  7  shows  the  influence  of  the  steam  to  carbon  ratio  on 
the  system  efficiency.  The  optimum  steam  to  carbon  ratio 
(1.8)  leads  to  a  maximum  efficiency  of  46%.  In  Fig.  7,  we 
have  located  the  self-sufficient  limit.  Here,  the  interest  of  the 
additional  firing  is  demonstrated:  without  supplementary 
firing,  the  maximum  efficiency  would  be  of  42%.  It  becomes 
46%  with  supplementary  firing.  The  maximum  performance 
for  this  operating  conditions  corresponds  to  an  additional 
firing  of  4%  based  on  the  total  fuel  inlet. 

3.2.  The  SMR  reactant  preheating  train 

From  the  analysis  of  the  composite  curves  of  the  system, 
several  improvements  might  be  suggested.  In  order  to  reduce 
the  plateau  at  the  steam  vaporization  level,  it  might  be 
envisaged  to  install  a  special  preheating  device  that  will 
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Composite  curves  of  the  fuel  cell  system  with  air  excess 


Fig.  6.  The  composite  curves  of  the  TPEMFC  system  with  a  steam  to  carbon  ratio  of  1.4  (SMR  and  PEMFC  temperatures  fixed  at  750  and  85  °C)  and  an 
oxygen  content  in  the  flue  gases  of  6%. 


saturate  the  natural  gas  with  a  liquid  water  injection  while 
preheating  the  mixture.  This  device  will  profit  from  the 
partial  pressure  effect  to  vaporize  part  of  the  water,  allowing 
to  transform  the  constant  temperature  plateau  into  a  tem¬ 
perature  profile  as  indicated  in  Fig.  8.  In  this  way,  the  cold 
and  hot  composites  are  closer  meaning  that  both  the  required 


energy  and  energy  losses  are  minimized.  Fig.  9  points  out 
the  comparison  between  the  efficiencies  with  and  without 
vaporization.  The  maximum  efficiency  of  the  system  is 
increased  of  about  3%  and  corresponds  to  a  higher  steam 
to  carbon  ratio.  The  efficiency  seems  to  stabilize  after 
the  maximum  is  reached  for  higher  steam  to  carbon  ratio. 


Fig.  7.  The  efficiency  as  a  function  of  steam  to  carbon  ratio  (SMR  and  PEMFC  temperatures  fixed  at  750  and  85  °C  and  a  fuel  utilization  of  95%). 
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Fig.  8.  The  composite  curves  of  the  TPEMFC  system  without  vaporization 
for  a  steam  to  carbon  ration  of  2.4  (SMR  and  PEMFC  temperatures  fixed  at 
750  and  85  °C  and  a  fuel  utilization  of  95%). 

The  limit  of  the  self-sufficient  process  is  pushed  to  the  right 
but  additional  firing  is  still  needed  to  reach  the  optimum. 
The  self-sufficient  system  has  a  higher  efficiency  when  the 
saturation  device  is  used. 

3.3.  The  air  combustion  and  depleted  fuel  preheating 

Air  combustion  and  depleted  fuel  preheating  (combustion 
preheating)  recovers  heat  from  below  the  reforming  tempera¬ 
ture  to  make  it  available  above  the  reforming  temperature. 
This  reduces  the  additional  fuel  required  and  increases  the 


self  sufficient  limit.  The  results  shown  in  Fig.  10  demonstrate 
that  the  optimum  efficiency  is  reached  for  smaller  steam  to 
carbon  ratio  (2.4  instead  of  4).  The  systems  with  preheating 
will  therefore  be  more  compact  and  less  expensive.  The 
maximum  efficiency  is  of  49%  with  preheating  while  it 
is  of  47%  without  preheating.  As  it  will  be  shown  in 
Fig.  15,  the  difference  in  efficiency  is  even  higher  for  lower 
fuel  utilization.  Working  with  combustion  preheating  will 
also  help  in  water  management  issues:  increasing  the  air 
excess  to  reduce  the  water  content  at  the  cathode  will  have 
a  smaller  effect  on  the  overall  performance  of  the  system 
(Fig.  11). 

3.4.  The  influence  of  the  steam  methane  reformer 
temperature  (T0j smr) 

The  developed  model  gives  information  about  the  optimal 
reforming  temperature.  In  Fig.  12,  we  present  the  efficiency 
as  a  function  of  the  steam  to  carbon  ratio  for  different  values 
of  To, SMR • 

Similar  optimal  system  efficiencies  may  be  reached  but  for 
different  steam  to  carbon  ratio.  For  similar  efficiencies,  the 
curve  corresponding  to  the  temperature  of  750  °C  corres¬ 
ponds  to  the  smallest  steam  to  carbon  ratio  and  to  the  highest 
efficiency.  7),  smr  =  750  °C  seems  to  be  the  optimal  tem¬ 
perature.  We  should  observe  that  the  curves  corresponding  to 
To, smr  =  700  and  800  °C  are  practically  supeiposed  mean¬ 
ing  that  a  range  of  T0  smr  will  provide  similar  efficiencies.  In 
Fig.  12,  we  also  compare  the  hydrogen  molar  fraction  at  the 
fuel  processing  subsystem  outlet  and  the  system  efficiency. 


—♦—With  vaporisation  —■—Without  vaporisation 


Fig.  9.  The  efficiency  with  and  without  vaporization  as  a  function  of  steam  to  carbon  ratio  (SMR  and  PEMFC  temperatures  fixed  at  750  and  85  °C  and  fuel 
utilization  95%). 
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—♦—without  preheating  -■—With  preheating 


Fig.  10.  The  efficiency  with  and  without  combustion  preheating  as  a  function  of  steam  to  carbon  ratio  (SMR  and  PEMFC  temperatures  fixed  at  750  and 
85  °C  and  a  fuel  utilization  of  95%). 


The  optimum  steam  to  carbon  ratio  for  the  fuel  processing 
and  for  the  global  efficiency  do  not  coincide  indicating  the 
importance  of  adopting  an  integrated  approach  rather  then 
separated  optimizations. 


3.5.  The  influence  of  the  cell  temperature  rceii 

Fig.  13  shows  the  relative  humidity  at  anode  and  cathode 
outlets  and  inlets  for  two  different  temperature  of  the  cell. 


— a — Anode  humidity  at  80% 
-  -♦ — Anode  humidity  at  95% 
— A —  Efficiency  at  80% 

♦  Efficiency  at  95% _ 


Fig.  11.  The  influence  of  the  fuel  utilization  on  the  self-sufficient  limit. 
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MF  at  Tref=800’G 
— ♦  - -MF  at  Tref=750°C 

-  -  -a-  -  MF  at  Tref=700°C 

-  -  -x  -  -  ■  MF  at  Tnef=650*C 
— • —  Eff  at  Tref=800“C 
— • —  Eff  at  Tref=750°C 
— * —  Eff  at  Tref=700°C 
— x —  Eff  at  Tref=650°C 


Fig.  12.  The  efficiency  and  the  H?  molar  fraction  for  different  reforming  temperatures,  the  cell  temperature  is  fixed  at  85  °C  and  a  fuel  utilization  of  95%. 


The  steam  to  carbon  ratio  and  the  reforming  temperature 
have  a  big  influence  on  the  humidity  and  therefore  on  the 
water  management  strategy.  As  all  results  have  been  carried 
out  with  the  mass  flow  rate  of  water  fixed  at  zero  through  the 
polymer  membrane,  a  careful  analysis  should  be  performed 
to  study  these  phenomena  and  introduce  the  conclusions  as 


constraints  in  the  simulation  and  integration  model.  The 
modeling  strategy  adopted  here  will  help  in  computing  the 
impact  at  the  system  level. 

The  cell  temperature  does  not  really  affect  the  overall 
system  efficiency.  The  impact  will  be  rather  on  the  tech¬ 
nology.  The  humidity  decreases  with  an  increase  of  the  cell 


— -  Cathode  output  at  75*0 
— • . . .  Cathode  output  at  85°C 
-  -  -  A-  -  Anode  output  at  75°C 
— x  -  -  Anode  output  at  85°C 
— ♦—  Cathode  input  at  75°C 
— • — Cathode  input  at  85°C 
— * — Anode  input  at  75°C 
— X — Anode  input  at  85°C 


Fig.  13.  Anode  and  cathode  humidities  at  the  inlet  and  outlet  of  the  cell  as  a  function  of  the  steam  to  carbon  ratio  for  two  different  cell  temperatures. 
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-  -  -x  -  -  Tref=800°C  at  95%  -  -  -x  -  -  Tref=750°C  at  95%  ------  ■Tref=700°C  at  95%  -  -  -a-  -  Tref=650°C  at  95% 

— x — Tref=800°C  at  80%  - x — Tref=750',C  at  60% - • - Tref=700°C  at  80%  — * — Tref=650X  at  80% 


Fig.  14.  The  efficiency  for  different  fuel  utilization’s  and  different  reforming  temperatures. 


temperature,  as  shown  in  Fig.  13,  reasonable  efficiency  can 
be  reached  within  the  optimal  operating  temperature  range 
and  with  acceptable  humidity.  The  difference  between  the 
inlet  and  the  outlet  humidity  is  higher  for  lower  temperatures 


(Fig.  13)  indicating  that  working  with  higher  temperatures 
will  create  less  problems  of  water  repartition  in  the  cell. 

According  to  the  affirmations  above,  the  cell  temperature 
has  to  be  as  high  as  possible  meaning  also  a  higher  energy 


-•—With  preheating  — Without  preheating 


Fig.  15.  The  efficiency  for  a  fuel  utilization  of  80%  with  and  without  preheating. 
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content  in  the  available  heat  for  combined  heat  and  power 
purposes. 

3.6.  The  influence  of  the  fuel  utilization 
in  the  cell  ( t]comh ) 

The  fuel  utilization  in  the  cell  is  an  important  factor  that 
affects  the  balance  of  the  plant  and  therefore  the  global 
efficiency  of  the  system.  A  higher  fuel  utilization  leads  to  a 
decrease  of  the  available  energy  for  post-combustion  leading 
therefore  to  the  possible  increase  of  the  additional  fuel 
requirement.  This  may  be  compensated  by  an  increase  in 
the  fuel  cell  power  production.  The  process  is  less  influenced 
by  the  fuel  utilization  with  the  higher  reforming  temperatures 
(Fig.  14).  The  difference  decreases  when  the  steam  to  carbon 
ratio  increases.  In  fact,  for  values  above  a  steam  to  carbon 
ratio  of  3,  the  curves  are  supeiposed.  For  lower  reforming 
temperatures,  the  situation  is  different  and  the  differences  are 
much  bigger. 

At  a  constant  steam  to  carbon  ratio  of  1.5,  the  difference 
reaches  up  to  10%  in  terms  of  system  efficiency.  For 
higher  steam  to  carbon  ratio,  the  difference  is  attenuated. 
The  importance  of  the  fuel  utilization  rate  is  higher  for 
^o,smr  =  750  °C  than  for  800  °C.  The  fuel  utilization  has 
also  an  influence  on  the  self-sufficient  limit  since  the 
depleted  fuel  has  a  lower  hydrogen  content  for  higher 
fuel  utilization  rate.  This  is  shown  in  Fig.  11.  A  lower  fuel 
utilization  does  not  need  additional  firing  to  reach  the 
maximum  efficiency  of  the  system,  but  in  this  case  the 
steam  to  carbon  ratio  is  much  higher.  This  might  simplify 
the  design  and  the  control  of  the  system.  One  should  remind 


that  these  calculations  are  made  with  a  stoechiometric 
combustion.  Increasing  the  air  excess  will  lower  the  self- 
sufficient  limit. 

The  influence  of  the  air  and  depleted  fuel  preheating  has 
been  analyzed  in  Fig.  15  for  a  fuel  utilization  rate  of  80%. 
This  figure  should  be  compared  with  Fig.  10  that  refers  to  a 
fuel  utilization  of  95%.  The  efficiency  gain  obtained  by 
preheating  the  depleted  fuel  and  the  combustion  air  to  a 
temperature  of  500  °C  is  of  9  points  of  efficiency  (from  39  to 
48%)  which  is  quite  considerable.  The  air  and  depleted  fuel 
preheating  heat  exchange  before  combustion  is  therefore  an 
important  device  in  order  to  maximize  the  overall  efficiency 
of  the  system  especially  when  lower  fuel  utilization  are 
expected  in  the  cell.  The  preheating  will  compensate  the 
efficiency  losses. 

4.  Conclusions 

A  combined  simulation  and  process  integration  model  has 
been  developed  to  optimize  the  performances  of  a  proton 
exchange  membrane  fuel  cell  system  for  stationary  applica¬ 
tions.  The  process  integration  techniques  have  been  used  to 
model  the  heat  exchanger  network  without  having  to  con¬ 
sider  its  configuration.  This  approach  reveals  to  be  very 
interesting  in  order  to  study  the  impact  of  the  major  deci¬ 
sions  parameters  (steam  to  carbon  ratio,  reforming  and  cell 
temperatures).  The  interest  of  the  process  integration  tech¬ 
niques  stands  in  the  fact  that  it  is  able  to  identify  heat 
exchange  bottlenecks  in  the  system  (pinch  points)  and  to 
propose  process  modifications  that  will  increase  the  overall 


Reforming  temperature 


— * — Fuel  utilization=95%  and  Tcell=85°C 
-  -  -X-  -  Fuel  utilization=95%  and  Tcell=75°C 
— a — Fuel  utilization=80%  and  Tcell=85°C 
— A-  -  Fuel  utilization=80%  and  Tcell=75°C 


Fig.  16.  Maximum  system  efficiency  as  a  function  of  the  reforming  temperatures  for  different  cell  temperatures  and  fuel  utilizations. 
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efficiency  of  the  system  (saturation/preheating  of  the  natural 
gas  feed,  combustion  preheating).  The  combined  use  of 
process  modeling  using  commercial  software  and  process 
integration  techniques  implemented  into  an  optimization 
frame  work  leads  to  propose  a  new  integrated  system 
structure  whose  optimum  is  not  the  optimum  of  the  sepa¬ 
rated  subsystem. 

The  model  has  been  used  to  optimize  the  system  perfor¬ 
mances  and  analyze  the  sensitivity  of  the  decision  para¬ 
meters  with  respect  to  the  major  unknowns  of  the  future  cell 
performance  (e.g.  fuel  utilization  in  the  cell).  The  results  are 
not  only  optimal  decision  parameters  values  but  also  optimal 
process  configurations.  In  Fig.  16,  we  present  the  evolution 
of  the  optimal  efficiency  as  a  function  of  the  steam  methane 
reformer  temperature  for  different  values  of  the  fuel  utiliza¬ 
tion  and  the  cell  temperature.  Each  point  in  the  curves 
correspond  to  an  optimal  value  of  the  steam  to  carbon  ratio 
and  to  a  given  process  configuration  (i.e.  an  investment)  that 
remains  to  be  determined  more  accurately  by  designing  the 
corresponding  heat  exchange  system  and  modeling  the 
integrated  system. 

The  combined  model  has  also  been  used  to  study  the 
influence  of  the  technological  constraints  to  be  considered 
for  the  cell  (e.g.  the  water  management  issues,  fuel  utiliza¬ 
tion)  and  deduce  the  corresponding  optimal  process  config¬ 
urations. 

The  solutions  obtained  from  the  application  of  the  method 
to  the  TPEMFC  system  feature  total  efficiency  as  high  as 
49%  that  are  much  higher  than  the  one  of  the  original  system 
design. 

The  results  of  this  preliminary  analysis  will  be  used  to 
design  the  heat  exchange  network  for  the  system  and  per¬ 


form  a  more  accurate  thermo-economic  optimization  where 
the  trade-off  between  the  efficiency  increase  w.r.t.  invest¬ 
ment  will  be  considered.  The  advantage  of  the  adopted 
approach  stands  in  the  fact  that  by  using  process  integration 
techniques  the  changes  in  the  heat  exchange  system  con¬ 
figuration  are  allowed  and  the  structure  will  adapt  according 
to  constraints  trade-off. 
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